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A self-consistent renormalization (RG) scheme has been applied to nonhelical magnetohydrody- 
namic turbulence with normalized cross helicity a c = and a c — > 1. Kolmogorov's 5/3 powerlaw is 
assumed in order to compute the renormalized parameters. It has been shown that the RG fixed 
point is stable for d > d c ~ 2.2. The renormalized viscosity v* and resistivity rf have been cal- 
culated, and they are found to be positive for all parameter regimes. For a c — and large Alfven 
ratio (ratio of kinetic and magnetic energies) ta, v* = 0.36 and rf = 0.85. As ta is decreased, 
v* increases and rf decreases, untill rA ~ 0.25 where both v* and rf are approximately zero. For 
large d, both v* and rf vary as d -1 / 2 . The renormalized parameters for the case a c — + 1 are also 
s ! . reported. 
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I. INTRODUCTION 



The Renormalization Group (RG) technique provided an important tool for understanding of the phenomena 
of phase transitions and critical phenomena. Motivated by this success, researchers have applied this technique 
to turbulence and other nonequilibrium systems. In this paper we apply RG to magnetohydrodynamic (MHD) 
^ \ turbulence. 

Forster et al. Q first applied dynamical RG procedure to the analysis of Navier-Stokes and Burgers equations, 
both stirred by random forces. They eliminated the large-wavenumber modes and included the effects of nonlinearity 
in the effective viscosity. Later, DeDominicis and Martin ||, Fournier and Frisch [3], Yakhot and Orszag Q, many 
others studied various aspects of RG for fluid turbulence, but all of them considered certain form of external forcing. 
McComb and his coworkers ( ||U and reference therein) instead applied self-consistent RG procedure; here the energy 
spectrum was assumed to be Kolmogorov's spectrum, and the renormalized viscosity was computed iteratively. We 
have adopted a similar scheme for MHD turbulence. 

Fournier et al. |?J , Camargo and Tasso |jj , and Liang and Diamond || employed RG technique to MHD turbulence 
on the similar lines as adopted by Forster et al. for fluid turbulence. Fournier et al. found that for d> d c ~ 2.8, 
there are two nontrivial regimes: a kinetic regime where the renormalization of the transport coefficients is due to 
kinetic small-scales; and a magnetic regime where it is due to the magnetic small-scales. In dimensions 2 < d < d c , 
there is no stable fixed point for turbulence with sufficiently strong external currents. Camargo and Tasso [pj studied 
the effective (renormalized) viscosity and resistivity and concluded that negative viscosity and resistivity are not 
permissible except in certain cases. Liang and Diamond || showed that no RG fixed point exists in 2D MHD 
turbulence. All the above authors however do not give any clear idea about the energy spectrum of MHD turbulence, 
which still remains controversial. However, Verma fnj used technique similar to that of McComb and Coworkers ( 
HU an d references therein) and showed that Kolmogrov's spectrum is a consistent solution of MHD RG equations. 

Kraichnan [0 and Irosnikov [Jl2| first gave phenomenology of steady-state, homogeneous, and isotropic MHD 
turbulence. They argued that the kinetic and magnetic energy spectra (E u (k) and E b (k) respectively) are 

E u (k) = E b (k) = A(nB ) 1/2 fc- 3/2 (1.1) 

where n is the total energy cascade rate, Bo is the mean magnetic field or the magnetic field of the largest eddy, and 
A is a universal constant of order one. Later Matthaeus and Zhou |13[] , and Zhou and Matthaeus proposed a 
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generalized phenomenology in which the cascade rates IT 11 and energy spectra £ ,± (fc) of Elsasser variables z ± (= u ± b) 
are related by 



, _ A 2 E+(k)E-(k)k 3 



IP = v w (1.2) 

B + y/kE±(k) 



where A is a constant. In the limit Bq = 0, we can immediately obtain 



(n=F) 2 / 3 



where _ftT ± are Kolmogorov's constants for MHD. The Eq. ( |l.3| ) was first derived by Marsch |l5[ . 

The energy spectrum of the solar wind, a well known observational platform for MHD turbulence, is closer to A: -5 / 3 
than to £;~ 3 / 2 [ ^6|JT7[| . In addition, the recent numerical 18-20) and theoretical []Tc|fjl|]2"2] | work support Kolmogorov- 



like phenomenology for MHD turbulence. Curiously, Kolmogorov's spectrum is observed even in the presence of a 



mean magnetic field. To understand theoretically, Verma |10|, Sridhar and Goldreich |2l|| , and Goldreich and Sridhar 
p2| applied field theoretic techniques to MHD turbulence. Verma |l(| showed that the mean magnetic field Bq 
responsible for the Alfv en e ffect gets renormalized in strong turbulence (B (k) oc A: -1 / 3 ). When we substitute the 



renormalized B in Eq. (1.1), the energy spectrum turns out to be proportional to fc -5 / 3 . Sridhar and Goldreich pl[ 
and Goldreich and Sridhar used wave interaction picture and showed that the energy spectrum is anisotropic, 
and it follows Kolmogorov's powerlaw for strong turbulence. 

In this paper we calculate the renormalized viscosity and resistivity for MHD. Here we test whether Kolmogorov's 
spectrum is a consistent solution of RG or not if the turbulence is forced at large-scale. We carry out wavenumber 
elimination and obtain recursive RG equations. Kolmogorov's spectrum is substituted for the correlation function 
in the recursive RG equation. After that we attempt to obtain a convergent solution for the renormalized viscosity 
and resistivity. If the convergent solution exists, then the Kolmogorov's spectrum is one of the correct choices for 
the energy spectrum of MHD turbulence. For these cases we compute many important quantities, e.g., renormalized 
viscosity, renormalized resistivity, cascade rates, etc., in a reasonably simple manner. In this paper we report the 
calculation of renormalized viscosity and resistivity, and in the subsequent paper (referred to as paper II) we will 
report the cascade rate calculation. 

There are two field variables in MHD: the velocity fluctuation u and the magnetic field fluctuation b. The quantities 
of interest are Kinetic energy (KE), magnetic energy (ME), total energy (KE+ME), and cross helicity (H c = u • b). 
The dimensionless parameters used in this paper are the ratio of twice the cross helicity and the total energy, called 
the normalized cross helicity a c (a c — 2H C /(KE + ME)), and the ratio of KE and ME, called the Alfven ratio ta- 
In this paper we have assumed that the mean magnetic field is zero, and so are magnetic and kinetic helicities, i.e., 
we are considering nonhelical plasma. The calculation of the renormalized parameters is quite complex for arbitrary 
<7 C and ta- Therefore, we limit ourselves to two limiting cases: (1) a c = and whole range of r^; (2) a c — > 1 and 
ta = 1- Even though full range of er c is observed in terrestrial and extraterrestrial plasmas, the case a c « is most 
pronounced. 

The outline of the paper is as follows: in sections 2 and 3 we compute the renormalized viscosity and resistivity for 
the cases a c = and a c — > 1 respectively. Section 4 contains the summary and conclusions. 

II. RG CALCULATION: a c = 

In this section we will calculate the renormalized viscosity and resistivity of MHD equations for zero mean magnetic 
field and zero cross helicity (a c = 0). In this case, the RG equations in terms of velocity and magnetic field variables 
are simpler as compared to those with Elsasser variables. Therefore, we will work with u and b variables. We take 
the following form of Kolmogorov's spectrum for KE [E u (k)] and ME [E b (k)] 

E u {k) = K u n 2/3 k~ 5/3 , (2.1) 
E b (k)=E u (k)/r A , (2.2) 

where K u is Kolmogorov's constant fo r MH D turbulence, and n is the total energy flux. In the limit a c = 0, we have 
E+ = E~ and n+ = n~ = n [cf. Eq. (pll- Therefore, E total (k) = E+(K) = E u (k) + E b (k) and 

K+ = ^"(l + r^ 1 ) (2.3) 
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With this preliminaries we start our RG calculation. The incompressible MHD equations in the Fourier space are 



-iu + vk 2 ) Ui(k) = -|Pi} m ( k ) / dp[u ] (p)u m (k -p)-b j {p)b m (k-p)} (2.4) 
(-iw + V k 2 ) k(k) = -iF£ m (k) / dp[u 3 {p)b m {k - p)] (2.5) 



fe«<(k) = (2.6) 

fe&i(k) = (2.7) 

with 

PfjmW = fejAm(k) + fc ro ^(k); (2.8) 

Pim(k) = S im - (2.9) 

^T m (k) = fc^ m - fc m %; (2.10) 

fc=(k,w); (2.11) 
dp = dpduj / {2n) d+1 . (2.12) 

Here 1/ and r\ are the viscosity and the resistivity respectively, and d is the space dimension. 

In our RG procedure the wavenumber range (ftjv,fco) is divided logarithmically into N shells. The nth shell is 
(k n , k n -\) where k n = h n ko(h < 1). In the following discussion, we carry out the elimination of the first shell (hi, ho) 
and obtain the modified MHD equations. We then proceed iteratively to eliminate higher shells and get a general 
expression for the modified MHD equations. The renormalization group procedure is as follows: 

1. We divide the spectral space into two parts: 1. the shell (fc l5 k ) — k y , which is to be eliminated; 2. (k N ,ki) = 
/c < , set of modes to be retained. Note that vq and t]q denote the viscosity and resistivity before the elimination 
of the first shell. 



2. We rewrite Eqs. (ETL 2.1) for k K and fc > . The equations for uf(k) and bf{k) modes are 



+ E^(fc) u<(k) + S^(fc)6f (fc) = -x^+Jk) / dp([uf(p)u<(k-p)] 



2 

+2{uf{p)u>(k-p)\ + [u>(p)u>(k-p)} 

—Similar terms for b) (2-13) 

(-iw + Egj(fc)) b<(k) + S^)(fc)u<(fc) = -iPy m {V) J dp([uf(p)b<(k - p)} 

+[uf(p)b>(k -p) + u>(p)b<(k - p)] 

+ [u>{p)b>(k-p)]) (2.14) 

The Es appearing in the equations are usually called the "self-energy" in Quantum field theory language. In 
the first iteration, E^t = v^k 2 and E^ = f]^k 2 , while the other two Es are zero. The equation for uf{k) 
modes can be obtained by interchanging < and > in the above equations. 



3. The terms given in the second and third brackets in the Right-hand side of Eqs. ( |2.13| , |2.14|) are calculated 
perturbatively. Since we are interested in the statistical properties of u and b fluctuations, we perform the usual 
ensemble average of the system j|. We assume that u > (fc) and b > (fc) have gaussian distributions with zero 
mean, while u<(fc) and h < (k) are unaffected by the averaging process. Hence, 



uf(k))=0 (2.15) 

6>(fc))=0 (2.16) 

uf(k))=uf(k) (2.17) 

b<(k)) =6<(A) (2.18) 
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and 

(u>(p)u>(q)) =P ij (p)C uu (p)S(p + q) (2.19) 

(b>(p)b>(q)) = P l0 {v)C bb {p)5{i> + q) (2-20) 

(u>(p)b>(q)) =P Z] (p)C mb (p)6(p + q) (2.21) 

The triple order correlations (^Xf {k)X> (p)X^{q)^j are zero due to Gaussian nature of the fluctuations. Here, X 

stands for u or b. In addition, we also neglect the contribution from the triple nonlinearity (^X<(k)X<(p)X<(q)^. 

The last assumption is used in many of the turbulence RG calculations Q|| . Refer to Zhou and Vahala pj| for 
a discussion on the importance of the triple nonlinearity. 

4. The details of the perturbative calculation is given in Appendix A. To the first order, the second bracketed 



term s of Eqs. (2.13, 2.14) vanish, but the nonvanishing third bracketed terms yield corrections to £s. Eqs. (2.13 



2.14) can now be approximated by 

(-iu> + E$ + £S$) u<(k) + (E$ + JS^) b<(k) = ~P$ m W J dp{u<(p)u<(k-p) 

-b<(p)b<(k-p)} (2.22) 

(-iw + E$ } + ££$>) &f (fe) + (E^ + uf (fc) = -iPy TO (k) | dp[uf(p)b<(k - p)} (2.23) 

with 

<5£|S(fc) = 73 L 2 / A #^(fc,p,g)G tl "(p)C«"(q) - 5 6 (fc,p, (? )G bb (p)C bb ( ( z) 

)p+q=k 

+S 6 (k lPl q)G ub (p)C ub (q) - S(fc,p,g)G b "(;3)G llb (<?)] (2.24) 
= (f) _\^ 2 [ Jp[-S{k,p,q)G uu {p)C ub {q) + S b {k, Pl q)G ub {p)C uu {q) 

+S(k lP ,q)G bu (p)C bb (q) - S 5 (k,p,q)G bb (p)C ub (q)} (2.25) 



(d 


- l)k 2 
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(d 


- l)k 2 




1 


(d 


- l)k 2 




1 


id 


- l)fc 2 



<5£^(fc) = ^— / _ #[5 8 (fc,p, (7 )G^(p)G tlb ( ( z) + S 10 {k,p,q)(^ b (p)C vb (q) 

+S 12 (k,p. / q)G ub (p)C bb (q) - S 7 (k,p,q)G bu (p)C uu (q)} (2.26) 
<5E bb } (fc) = (f} J Uh2 f Jp[-S s (k lP ,q)G uu (p)C bb (q) + S 9 (k lPl q)G bb (p)C uu (q) 

" "" / p+q=k 

+S 11 (k,p,q)G ub (p)C ub (q) - S 9 (k,p,q)G bu (p)C ub (q)} (2.27) 

The quantities Si(k,p, q) are given in the Appendix ||. The integral A is to be done over the first shell. 

The full-hedge calculation of Es is quite involved. Therefore, we take two special cases: (1) C ub — or a c = 0; 
and (2) C ub « C uu « C ub or a c — > 1. In this section we will discuss only the case a c = 0. The other case will 
be taken up in the next section. A word of caution is in order here. In our calculation the parameters used are 
a c (k) = 2C ub (k)/(C uu (k)+C bb (k)) and r A (k) = E u (k)/E b (k). These parameters differ from the global a c and 
ta, yet we have restricted ourselves to the limiting cases of constant a c (k) and rvt(A;) because of simplicity. 

When er c = 0, an inspection of the self-energy diagrams shows that E ub = £ b " = 0, and G ub = G bu = 0. Clearly, 
the equations become much simpler because of the diagonal nature of matrices G and E. The two quantities of 
]" " and 5E bb } 



interest S'Efff. and ST,^ will be given by 



A 



5£™(fc) = — - / dp (S(k,p,q)G uu (p)C uu (q) - S 6 (k,p,q)G bb (p)C bb (q)) (2.28) 

& J- Jp-\-q—k 

<5£ bb (fc) = -i- [ A dp(-Ss(k,p,q)G uu (p)C bb (q) + S 9 (k,p iq )G bb (p)C uu (q)) (2.29) 

1 Jp+q=k 
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The expressions for <5E's involves correlation functions and Green's function, which are themselves functions of 
S's. In earlier RG calculations (see e.g. the correlation function C{k) is function of forcing noise spectrum. 

In our self-consistent scheme we assume that we are in the inertial range, and the energy spectrum is proportional 
to Kolmogorov's 5/3 powerlaw. After that the renormalized Green's function, or renormalized v is computed 
iteratively. 

The frequency dependence of the correlation function are taken as: C uu (k,uj) — 2C UU (k)$l(G uu (k , u>)) and 
C bb (k,uj) — 2C bb (k)di(G bb (k,uj)). In other words, the relaxation time-scale of correlation function is assumed 
to be the same as that of corresponding Green's function. Since we are interested in large time-scale behaviour 
of turbulence, we take the limit u going to zero. Under these assumptions, the above equations become 

(d-l)k 2 J p+q=k (2tt)' 1 i/ (0 ) {v)V + ^(o) (<?)<? 

C bb (q) 

-S 6 (k,p,q) r — -g] (2.30) 

v(o) (p)p + V(o) im 



1 f d P r o _,u _ ^ C (l) 

/p+q= 



(d - l)k 2 J p+q=k (27r) d ^(o) (p)p^ + 77 (0) {q)q 2 



+S 9 (k,p,q) W 2 ] (2.31) 

v(o){p)p + v {o)(m 

Note that i/(Jfc) = T, uu (k)/k 2 and ^(fe) = H bb (k)/k 2 . 

There are some important points to remember in the above step. The frequency integral in the above is done 
using contour integral. It can be shown that the integrals are nonzero only when both the components appearing 



the denominator are of the same sign. For example, first term of Eq. (2.31) is nonzero only when both V(o)(p) 
and ?7(o) (q) are of the same sign. 

Let us denote v^{k) and ?7( 1 )(fc) as the renormalized viscosity and resistivity respectively after the first step of 
wavenumber elimination. Hence, 

"(i)(*)=f(0)(*)+^(0)(*); (2-32) 
V(i)(k)=V(o)(k)+5 V(0) (k) (2.33) 

5. We keep eliminating the shells one after the other by the above procedure. After n + 1 iterations we obtain 

f(n+i) ( fc ) = v (n) {k) + 5v( n ) (k) (2.34) 
V(n+i) (k) = T)(n) (k) + 5i] (n) (k) (2.35) 



where the equations for 5u/ n -\(k) and 5rj^{k) are the same as the Eqs. (2.30, 2.31) except that V(n\(k) and 
?7(o)(fc) appearing in the equations are to be replaced by i/( n )(k) and r)r n -)(k) respectively. Clearly vi n +\){k) and 
T /(n+i)(k) are the renormalized viscosity and resistivity after the elimination of the (n + l)th shell. 

6. We need to compute 6u^ and 8r)/ n \ for various n. T hese comp utations, however, require vr n \ and T)/ n y In our 
scheme we solve these equations iteratively. In Eqs. ( [2.30 , 2.31) we substitute C(k) by one dimensional energy 
spectrum E(k) 

C uu^ = 2(2tt) k -(d-D E u^ ( 2 36) 

C bb {k) = s 2 ^. ' k~^E b (k) (2.37) 



where Sd is the surface area of d dimensional spheres. We assume that E u (k) and E b (k) follow Eqs. (O O) 
respectively. Regarding vr n -\ and r}r n -\ , we attempt the following form of solution 

u (n) (k n k>) = (K^U^k-^^ik') (2.38) 

V(n) (k n k>) = (in 1/2 n 1/3 fc- 4/3 ^)(fc') (2-39) 
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with k = k n+ ik'(k' < 1). We expect v^(k') and rfL\{k') to be universal functions for large n. The substitution 
of C uu (k), C bb (k), vi n \{k), and rjt n \{k) yields the following equations: 



1 



S(k',p', q 



, 2 E u (g>) 



1 



S 6 (k', 



r 7( *„ ) (^'y 2 + ^„)(V)9' 2 



(rf-l)y p ' +q '= k '^ (d-l)5«, g" 1 - 1 

1 



-S 8 (fc',pW)- 



^ n) (hp')p' 2 + r,* in) (hq')q' 2 
v* n+ i)(k') = h 4 / 3 u* {n) (hk') + ft- 4 / 3 ^ ( * n) (fc') 
^ + i)(fc') = h 4/3 r)* n) ( hk ') + h- 4/3 6 v l n) (hk') 



+ S Q (k',p',q') — 



r,t n) (hp')p' 2 +v* {n) (hq')q' 2 



(2.40) 



(2.41) 

(2.42) 
(2.43) 



where the integrals in the above equations are performed iteratively over a region 1 < p',q' < 1/h with the 
constraint that p' + q' = k'. Fournier and Frisch || showed the above volume integral in d dimension to be 



f dq! = f dpdq (^Y 2 (sma) d - 

Jp'+q'=k' J \ K / 



(2.44) 



where a is the angle between vectors p and q. 

7. Now we solve the above four equations self consistently for various ta$- We have taken h = 0.7. We start with 
constant values of v*^ and 77^ and compute the integrals using Gauss quadrature technique. Once Sv* ^ and 
<5r7* ^ have been computed, we can calculate i/JU and 77?^. We iterate this process till ^(* m+1 )(fc') ~ v * m ){k') 
and V( m +i)(k') ~ ^(m)^')) that is, till they converge. We have reported the limiting v* and 77* whenever the 
solution converges. The criterion for convergence is that the error must be less than 1%. This criterion is usually 
achieved by n = 10 or so. The result of our RG analysis is given below. 

We have carried out the RG analysis for various space dimensions and find that the solution converges for all 
d > d c 2.2. Hence, the RG fixed-point for MHD turbulence is stable for d > d c . For illustration of convergent 
solution, see Fig. [l] for the plots of v* n ^(k') and rj^(k') for d = 3, ta = 1- The RG fixed point for d < d c is unstable. 
Refer to the plots of Fig. |2| for d = 2, ta = 1 as an example of an unstable solution. From this observation we can 
claim that Kolmogorov's powerlaw is a consistent solution of MHD RG equations at least for d > d c . The values of 
asymptotic (k' — > limit) v* and 77* for various d and Ta = 1 are displayed in Table [j| Clearly v* and 7/* decreases 
with the increase in d for d > 4. The plot of Fig. (||) shows that both v* and 77* are proportional to dr 1 ! 2 . In the 
same plot we have also displayed v* for pure fluid turbulence; there too v* oc d~ x / 2 . This result is consistent with 
Fournier et al.'s prediction for fluid turbulence |23]| . 

An important point to note is that the stability of RG fixed point in a given space dimension depends on Alfven 
ratio and normalized cross helicity. For example, for d = 2.2 the RG fixed point is stable for ta > 1, but unstable for 
ta < I- A detailed study of stability of the RG fixed point is required to ascertain the boundary of stability. 

The values of renormalized parameters for d — 3 and various ta are shown in Table ||. For large ta (fluid dominated 
regime), v* is close to renormalized viscosity of fluid turbulence (ta — 00), but 77* is also finite. That means there is 
a positive magnetic energy flux even when ta — > 00. As ta is decreased, 77* decreases but v* increases, or the Prandtl 
number Pr = vjr\ increases. This trend is seen till ta ~ 0.25 where the RG fixed point with nonzero v* and 77* 
becomes unstable, and the trivial RG fixed point with v* = 77* = becomes stable. This result suggests an absence 
of turbulence for ta below 0.25 (approximately). Note that in the ta — * (fully magnetic) limit, the MHD equations 
become linear, hence there is no turbulence. Surprisingly, our RG calculation suggests that turbulence disappear near 
r A = 0.25 itself. 

The final v*{k') and T]*(k') are constant for small k' but shifts towards zero for larger k' (see Fig. |l|). Similar 
behaviour has been seen by McComb and coworkers |p3 for fluid turbulence; this behaviour is attributed to the neglect 
of triple nonlinearity, and the corrective procedure has been prescribed by Zhou and Vahala p5| . For simplicity we 
have not included the effects of triple nonlinearity in our calculation. 
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Kraichnan's k 3 ' 2 energy spectrum [see Eq. ([H])] and a uu = a bb oc kBo do not satisfy the RG equations. This 
implies that Kraichnan's 3/2 powerlaw is not a consistent solution of RG equations. 

Pouquet [^6| and Ishizawa and Hattori |27]] calculated the contributions to the renormalized (eddy) viscosity and 
resistivity from each of the four nonlinear terms of MHD equations using EDQNM (Eddy-Damped-Quasi-Normal- 
Markovian) approximation. This calculation has been done for d — 2. Let us denote the quantities u uu , v ub , rj bu , rj bb 
as contributions from the terms u • Vu, — b • Vb, — u • Vb, b • Vu respectively. Note that various cascade rates of 
MHD turbulence Q can be inferred from these parameters. The cascade rate from inside of the X sphere (X <) 
to outside of the Y sphere (Y >) is 

= / 2v XY {k)k 2 E x {k)dk (2.45) 
Jo 

where X, Y denote u or b. We have computed v XY for d> d c ; their values are listed in Tables | and ||. Our results 
show that r] bb is positive. Pouquet argues that rj bb is negative, while Ishizawa finds it to be positive. Since our RG 
procedure is applicable only for d > d c = 2.2, it is not proper to compare our results with those of Pouquet p(| and 
Ishizawa and Hattori p7ft . However, we can claim that the magnetic energy cascade rate (n£>) is positive for all 
d > d c because i] bb > 0. 

The eddy viscosity and resistivity listed in Table | shows that j] bb * and v ub * decrease as d increases. On the contrary, 
rj bu * and v uu * first increases till around d = 5 then decreases as d increases. Since MHD fluxes are proportional to these 
parameters, the corresponding fluxes decrease as d is increased beyond 5. Near d — 2, v uu * is negative, reminiscent 
of inverse cascade energy in two-dimensional fluid turbulence in the 5/3 region. Also, rj bu * is negative near d = 2, i.e., 
there is a negative energy cascade from 6-sphere to outside u-sphere The above results are consistent with the recent 
simulation results of Dar et al. [^9| on two-dimensional MHD turbulence. 

The variation of v XY * with ta for d = 3 is displayed in Tables [n]. For large Ta, v wu * is close to the fluid turbulence 
v* . However, rj bb * also finite. Hence there is a significant magnetic energy cascade from 6-sphere to outside 6-sphere 
even in fluid dominated case. As ta is decreased, v ub * and rj bu * increase, hence cascade rates from u-sphere to outside 
6-sphere, and 6-sphere to outside u-sphere start becoming important. Incidently, both v uu * and rj bb * decrease as ta 
decreases. 

In this section we have calculated the renormalized viscosity and resistivity for a c = 0. In the next section we will 
calculate these parameters for the a c — > 1 limit. 



III. RG CALCULATION: a c -» 1 

The u and b correlation is very high when a c —*l. For this case it is best to work with Elsasser variables z ± = u ± b. 
In terms of Elsasser variables, ^|z + | 2 ) ^> ^z~| 2 ). These types of fluctuations have been observed in the solar wind 
near the Sun. However, by the time the solar wind approaches Earth, the normalized cross helicity is close to zero. 

In this section we will briefly discuss the RG treatment for the above case. For the following discussion we will 
denote (\z-\ 2 ) / (\z+\ 2 ) = r = (1 - a c )/(l + a c ). Clearly r < 1. 

MHD equations in terms of Elsasser variables are 

(-tw + v { o)±±k 2 ) zf{k) + m±T k 2 zf(k) = -iM ijm (k) J dp[zf(p)zt(k - p)] (3.1) 

where 

M ljm (k) = kjPimik) (3.2) 

Note that the above equations contain four dissipative coefficients v±± and v>± zf: instead of usual two constants 
v± = {y ± rf)/2. This is because -I — symmetry is broken in this case, consequently RG generates the other two 
constants. We carry out the same procedure as outlined in the previous section. After n + 1 steps of the RG 
calculation, the above equations become 

[-iu T (i>( n )±±(k) + 5v [n)±± (k)) k 2 ] zf<(k) 

+ (v(n)± T (k) + 5v (n)±T (k)) k 2 zj<{k) = -iM ijm (k) J dpzj < (p)zt l <{k^p) (3.3) 

with 
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A 



&>(„)++(*) = (a Tux / dp[S 1 (k,p, q )Gp ) (p)C--(q)+S 2 (k,p,q)G+T(p)C--(q) 

(a — L)K Jp + q =k 



+S 3 (k,p, q)G^ (p)C + - (q) + S 4 (k,p, q)Gfc (p)C + - (<?)] (3.4) 
5v {n)+ _ (k) = {d _\ )k2 _ k dp[S x (k,p, q)G+ n) (p)C- + (q) + S 2 (k,p, q)Gp ) (p)C- + (q) 

+Si{k 1 p,q)GJ- ) {p)C ++ {q) + S i (k,p 1 q)G- ( +{p)C ++ (q)] (3.5) 



where the integral is performed over the (n+l)th shell (k n+ i, k n ). The equations for the other two Siss can be obtained 

by interchanging + and — signs. Now we assume that the Alfven ratio is one, i.e., C H = E u — E b = 0. Under this 
condition, the above equations reduce to 

Ko++(*) = {d _\ )k 2 jT^_.rfp[5i(fc,P,g)gf n t(p) + g2(fc,P,g)Gf n 7(p)]g--(g) (3.6) 

Ho+-(*0 = {d _\ )k2 £ dp[S 3 (k,p, q)G^- ) ip) + S 4 (k,p, q)Gl+(p)]C ++ (q) (3.7) 

8v(n)- + (k) = {d \ )k2 ^_,^[^3(fc,P,g)Gf4(p) + ^4(fc,P,g)G^(p)]C--(g) (3.8) 

6v (n) —(k)= ( ,\ ]k 2 / JP[Si(k,p,q)G^(p) + S 2 (k,p,q)Gp ) (p)]C+ + (q) (3.9) 



The inspection of Eqs. ( p.6| -p.S|) reveal that v ++ and v |_ are of the order of r. Hence, we take the v matrix to be of 

the form 

*(*,«)= (3 ;) (3.io) 



It is convenient to transform the frequency integrals in Eqs. [(3.6-3.9)] into temporal integrals, which leads 

A dp 



5v [n)++ (k)= (d _ 1)fc2 J + k -^p J_ dt / [S 1 (k,p,q)G++(p,t-t') + S 2 (k,p,q)G+- ) (p,t-t')} 

C—(q,t-t') (3.11) 

and similar forms for equations for other vs. Green's function G(k, t — t') = exp — [i)k 2 (t — t')] can be easily evaluated 
by diagonalizing the matrix v. The final form of G(k, t — t') to leading order in r is 

^(1 - e-W-*')) - {f + f (| - (1 - e-/»(*-*')) 



G(M-f)= 7 , » 1/3 / V P JJ y , ' (3.12) 



The correlation matrix C(k,t — t') is given by 

fC++(k,t-t') C+-(k,t-t')\ , ( ■c++{k) C+-(k)\ 

yc-+(k,t-t') C—(k,t-t') ) - U ^> T r ^c*-+(fc) C— (k) J [ > 

The substitution of correlation functions and Green's functions yield the following expressions for the elements of 8v 

= t 7%^( q ){ Sl (k,p,q) 1 



(d - l)fc 2 J (2^) d WL 1V '"'^(g^ 2 

at( n )(p) / 1 1 



+S 2 (k,p,q)- 



0(n) (p) V 0{n) {P)P 2 + P(n) (g)g 2 P(n) (g)g 

~ S3{k,Piq) iM (p {n) (p)p 2 +P(n)(q)q 2 I3 {n) (p)p 2 ) } (3 ' 14) 

fo( " )(fc) = (^W/ wr S3iKp ' q) i^^ (3 - 15) 
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"(«)(<?) / 1 



+S 2 (k,p, q)- 



f3(n) (q) V P(n) (P)P 2 + P(n) (<l)<l 2 P(n) (p)P 



"(n)(p) 1_ 



+S 4 (k,p,q) a (n Z{ a ( , 2 ] (3-16) 



^"(Rp/ 4 ^^ (3 - 17) 

Note that 5a, <5/3, <5C, an d hence a, /J, £, tp, are all independent of r. To solve the above equations we substitute 
the following one-dimensional energy spectra in the above equations: 

E+(k) = K+± ^fc _5/3 (3.18) 

(n-) 2/3 

E~(k) =rE + (k), (3.19) 

For the elements of v we substitute 

Z (B) (Jfe) = ^ B) V^^-i i7S *-V3 (3.20) 

where Z stands for (,a,ip,(3. The renormalized Z*s are calculated using the procedure outlined in the previous 
section. For large n their values for d = 3 are 

% _ /0.86r 0.14 \ 

Z _ l^0.16r 0.84 J ' ^1) 

and for d — 2 they are 

~ ! I. ] Or ().") I 

Note that the solution converges for both d = 2 and <i = 3. 

As discussed in the earlier section, the cascade rates II ± can be calculated from the renormalized parameters 
discussed above. Using the energy equations we can easily derive the equations for the cascade rates, which are 

/■few 

n+ 



= / 2</c 2 ^ + (fc)+ / 2ak 2 {E u (k)-E h (k)) (3.23) 
Jo Jo 

/■few rkN 

n" = / 2(3k 2 E-(k)+ 2r^k 2 (E u (k)-E b (k)) (3.24) 
Jo Jo 



Under the assumption that rvt = 1, the parts of IT 1 proportional to (E u (k) — E b {k)) vanish. Hence, the total cascade 
rate will be 

n = i(n+ + rr) (3.25) 

fcjv 

(C + f3)k 2 E + {k) (3.26) 

Since £ and /3 are independent of r, the total cascade rate is proportional to r (for r small). Clearly the cascade rate II 
vanishes when r — or a c = 1 . This result is consistent with the fact that the nonlinear interactions are absent when 
only pure Alfven waves (z + or z") are present. The detailed calculation of the cascade rates II ± and the constants 
K ± is left for future studies (paper II). 
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IV. SUMMARY AND CONCLUSIONS 



In this paper we have constructed a self-consistent RG scheme for MHD turbulence and computed the renormalized 
viscosity and resistivity. These quantities find applications in simulations specially large-eddy-simulations (LES). 
In our RG scheme we assume that we are in the fully nonlinear range, and the renormalized energy spectrum is 
substituted for the correlation function. After the substitution of the correlation function, the renormalized Greens 
function or renormalized viscosity and resistivity are computed iteratively. In our procedure Kolmogorov's powerlaw 
is taken for the energy spectrum, and it is shown to be a self-consistent solution of the RG equation for d > d c sa 2.2. 
This result is consistent with the recent theoretical |lOU2l| , |22|| , numerical Jl^-]2C[] , and observational studies of solar 
wind that favor Kolmogorov's spectrum for MHD turbulence. Note that we do not carry out vertex (the 

coefficient of the nonlinear term) renormalization; there is an implicit assumption that the vertex renormalization, if 
any, is included when we substitute renormalized energy spectrum (Kolmogorov's) in the RG equation. 

For simplicity of the calculation we have taken two special cases: (1) a c = and full range of ta', (2) a c — > 1 
and Ta = 1. For these two cases, the renormalized viscosity and resistivity have been calculated for various space 
dimensions. For a c — 0, there exists a stable RG fixed point for d > d c w 2.2. The RG fixed point is unstable for 
d < d c . Our result is consistent with that Liang and Diamond (9j where they conclude that the RG fixed point for 
d = 2 is unstable. Note however that the stability depends on the value of ta and o~ c . An exhaustive study is required 
to ascertain the boundary of stability as a function of d, er c , and ta- 

Some interesting trends emerge as we vary space dimensionality. For ta = 1> the parameters v ub * and n bb * decrease 
with the increase of d, but v uu * and n bu * first increase till d ~ 5 then decrease. 

For d = 3, variation of ta shows some interesting features. For large ta, v* is close to the v* for pure fluid 
turbulence, but rf is also finite for this case. As ta is decreased, v* increases and n* decreases, or Pr — vjr\ increases. 
Another important result of our calculation is absence of turbulence for all ta below 0.25 or so. In the ta — > 
(fully magnetic) limit, the MHD equations become linear, hence there is no turbulence. It is however surprising that 
turbulence disappears near r a — 0.25 itself. 

Our RG calculation for fluid turbulence (?a — oo) show several interesting features. The RG fixed point is stable 
for d — 2, and the renormalized viscosity is approximately —0.5, a negative number. This is consistent with the inverse 
cascade of energy for the 5/3 region. Note however that for d — 3, the renormalized fluid viscosity is positive (0.38), 
hence the energy cascade is forward, consistent with Kolmogorov's hypothesis. 

The values of the renormalized viscosity and resistivity calculated by our calculation differs from that of Verma and 
Bhattacharjee j^o), the later calculation being similar to Kraichnan's Direct-Interaction- Approximation calculation 
for fluid turbulence. Verma and Bhattacharjee ]30|] assumed v ++ = v h and = v , which is correct only in a 
limited region of parameter space. In addition, they assumed a cutoff for the self-energy calculation to cure infrared 
divergence. In the present calculation we have overcome both these defects, and expect that the numbers reported 
here will match with the future simulation results. It is unfortunate that the RG scheme fails for d = 2, hence, several 
numerical results of 2D-MHD turbulence p^- |2p| , |29| could not be compared with the analytic results presented here. 

Kr aichnan's 3/2 energy spectrum and = E oc (kBo) do not satisfy renormalization group equations [Eqs. ( 2.30 - 
2.33)]. Hence E(k) oc £;~ 3 / 2 is not a consistent solution of RG equations. This result is in agreement with the recent 
calculation by Verma ]l0| , where it was shown that Bq = constant does not satisfy RG equations, but the renormalized 
mean magnetic field (Bq oc fc -1 / 3 ) and Kolmogorov's energy spectrum (/c~ 5 / 3 ) do satisfy the RG equations. From 
these arguments we can claim that Kraichnan's 3/2 powerlaw is ruled out for strong MHD turbulence. However, 
£-3/2 ma y s ^-jj k e considered in the framework of weak turbulence theory. 

In our calculation we take r^(fc) = E u (k)/E b (k) to be a constant. The global Alfven ratio E u /E b may differ 
significantly from r^(fc), still we make the approximation rA(k) = ta to simplify the calculation. Similarly we have 
assumed that a c (k) = const — a c . In addition we also assume isotropy which does not hold in the large ta limit. 
It is hoped that future calculations will be able to relax these assumptions and attempt solutions for more realistic 
situations. 

In many of the earlier RG calculations jjUl it is assumed that the dynamical system under consideration is forced 
by random noise at all scales. In those calculations, the noise and the vertex (coefficient of the nonlinear term) 
usually get renormalized along with the renormalization of dissipative coefficients. In contrast, in our scheme only 
the dissipative constants are renormalized, and the renormalization correction to other parameters (e.g. noise) are 
implicitly lumped into the full correlation function by assuming Kolmogorov's powerlaw. Our procedure is relatively 
simple, and one can easily calculate various interesting quantities. We hope that future developments in dynamical 
RG will be able to justify some of the assumptions made in our procedure. 

To conclude, we have applied RG procedure to MHD turbulence and calculated the renormalized viscosity and 
resistivity. These parameters find applications in simulations. We have shown that that Kolmogorov's spectrum is 
a consistent solution of MHD RG equation. The results presented here are in general agreement with the recent 
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EDQNM and numerical results. Validation of the results presented here with numerical simulations will give us 
important insights into physics of MHD turbulence. 
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APPENDIX A: PERTURB ATI VE CALCULATION OF MHD EQUATIONS 



The MHD equations (2.4, |2.5|) can be written as 



Ui(k) 
h(k) 



G uu (k) G ub {k) 
G bu (k) G bb {k) 



i^jm( k ) / dp[u 3 {p)u m (k - p) -bj{p)b m {k - §)} 
- iP ijm( k ) I dp[uj(p)b m (k - p)} 



(Al) 



where the Greens function G can be obtained from G 1 (k) 

G [k,u) = I s&u 



-iu - ~Z bb 



(A2) 



We solve the above equation perturbatively keeping the terms to the first order (nonvanishing). As usual, we 
represent the integrals by Feynmann diagrams. To the leading order, the quantities ttj and bi are expanded as 




Qi^b \b 



(A3) 




Quo \0„ 




G 



(A4) 



The variables u and b are represented by solid and dashed line respectively. The quantity G uu , G ub , G bb , G bu are 
represented by thick wiggly (photon), thin wiggly, thick curly (gluon), and thin curly lines respectively. The filled 
circle represents (—i/2)P^ m vertex, while the empty circle represents —iP^ m vertex. 



The substitution of uf and bf in the second bracketed terms of Eqs. (2.13, 2.14), denoted by IV, and l\ respectively, 
will yield 



T" 



r — 
1 2 — 




2< 



< 



(A5) 



(A6) 



We illustrate the expansion of one of the above diagrams: 



= 2 



0-+ 0-+ 2 €> 



<> 



<> 



<< 



« 
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<> 



<> 



<< 



<< 



rj > <^+ rj>>"+ r><r 

<> <> << << 



<> 



<> 



+higer oder diagrams 



<< 



« 



(A7) 




In the above diagrams solid lines denote C uu , and the dotted lines denote C . The correlation function C is 
denoted by dashed line. As mentioned earlier, the wiggly and curly lines deno te va ri ous G reen's functions. All the 
diagrams except 4,8,12, and 16th can be shown to be trivially zero using Eqs. ( 2.15| — 2.21 ). We assume that 4,8,12, 
and 16th diagrams are also zero, as usually done in turbulence RG calculations [Q]5| . Similarly we can show that all 
other diagrams of 1% and l\ are zero , henc e, 1^ = l\ = to first order. 

The third bracketed terms of ( 2.13 , 2.14 ), denoted by 1% and if respectively, are diagrammatically represented as 

- <: 



= -6E uu (k)- 



< 



■-6T, u \k)- 



< 



(A8) 



1 :>, 



O 



<5£ btl (fc)- 



< 



< 



(A9) 



where 



(d- 1)5^ 



-(d-l)<5£ b6 = 2 



>> 

>> 



-(d- l)<5£ ub = - 4 \^Jf 



-(d-l)SE bu = 2 1 



>> 



» 



» 




» 



2 

>> 
O Q 

» 




» 

2 *L P 

» 

- 4« • 

» 



» 



» 

» 
>> 

y 

>> 

6 Q 



(A10) 



(All) 



(A12) 



(A13) 



In Eqs. (A10-A13) we have omitted all the vanishing diagrams (similar to those appearing in J 2 ). Looking at 
the all the terms, we observe 1% contributes to £ ,m and £ ub renormalization, whi le Jo contr ibutes E bb and £ bu 
renormalization. The algebraic expressions for the above diagrams are given by Eqs. |^). 
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APPENDIX B: VALUES OF Si 



Si(k,p,q) are formed by contracting tensors PA m ,Pij m ,Pij and Mij m (see Eqs. (2.8,2.10 3~2])). 
SiS become functions of k,p, q and cosines x, y, z defined as 



p q = -pqx; 

The algebraic expressions for various S 1 ., 

Si(k,p,q 

S 2 {k,p,q 

S3(k,p,q 

Si{k,p,q 

S(k,p,q 

Sb{k,p,q 
S 6 (k,p,q 
S 7 (k,p,q 
S 8 (k,p,q 
Sg(k,p,q 
Sio(k,p,q 
S n (k,p,q 
S 12 (k,p,q 



There are many useful relationships between Si(k,p, q)s. Some of them are 



S(k,p,q) = S 
S5(k,p,q) = S 



q • k = qky; 



k,p, q) are given below. 



p • k = pkz. 



= M bjm (k)M rnab (p)P ja (q) 
= kp(d-2 + z 2 )(z + xy) 
= Maj m (k)M mab (p)P jb (q) 
= kp(—z + z 3 + y 2 z + xyz 2 ) 
= M bjm (k)Mj ab (p)P ma (q) 
= kp(—z + z 3 + x 2 z + xyz 2 ) 
= M ajm (k)M jab (p)P mb (q) 
= kp(—z + z 3 + xy + x 2 z + y 2 x + xyz 2 ) 

= p b + im (k)PL b (p)PM 

= kp ((d - 3)z + 2z 3 + (d- l)xy) 
= P b + jm (k)P mab (p)PM 
= kp ((d - l)z + (d- 3)xy - 2y 2 z) 

= Pa jm (k)P mba (p)P 3 b(q) 
= ~S 5 (k,p, q) 

= P^ m (k)PLb(p)PMPib(k) 

= S 5 (p,k,q) 

= P^m(k)P^ ab (p)Pma(q)Pi b (k) 

= -S 5 (p,k,q) 

= Pi~ jm (k)P- ab (p)Pja(q)P ib (k) 
= kp(d — l)(z + xy) 

= Pr xm {k)P mab {p)PAq)Pa{k) 

= S 9 (k,p, q) 

= P i 7 jm (k)P jab (p)P m a(q)P b (k) 
= -Sg(k,p, q) 

= ^ m (k)Pf ab (p)Pmb(q)Pia(k) 
= Sg(k,p, q) 



k,p, q) + S 2 (k,p, q) + S 3 (k,p, q) + 5 4 (fc,p, q) 
k,P, q) ~ S 2 (k,p, q) + S 3 (k,p, q) - S A (k,p,q) 



On simplification 
(Bl) 

(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
(B8) 
(B9) 
(BIO) 
(Bll) 
(B12) 
(B13) 
(B14) 



(B15) 
(B16) 



[1] D. Forster, D. R. Nelson, and M. J. Stephen, Phys. Rev. A 16, 732 (1977). 
[2] C. DeDominicis and P. C. Martin, Phys. Rev. A 19, 419 (1979). 
[3] J. D. Fournier and U. Frisch, Phys. Rev. A 17, 747 (1979). 



13 



[4] V. Yakhot and S. A. Orszag, J. Sci. Comput. 1, 3 (1986). 

[5] W. D. McComb, The Physics of Fluid Turbulence (Claredon, Oxford, 1990). 

[6] W. D. McComb, Rep. Prog. Phys. 58, 1117 (1995). 

[7] J. D. Fournier, P.-L. Sulem, and A. Pouquet, J. Phys. A 15, 1393 (1982). 

[8] S. J. Camargo and H. Tasso, Phys. Fluids B 4, 1199 (1992). 

[9] W. Z. Liang and P. H. Diamond, Phys. Fluids B 5, 63 (1993). 
[10] M. K. Vcrma, Phys. Plasma 6, 1455 (1999). 
[11] R. H. Kraichnan, Phys. Fluids 8, 1385 (1965). 
[12] P. S. Iroshnikov, Sov. Astron. I. 7, 566 (1964). 
[13] W. H. Matthaeus and Y. Zhou, Phys. Fluids B 1, 1929 (1989). 
[14] Y. Zhou and W. H. Matthaeus, J. Geophys. Res. 95, 10291 (1990). 

[15] E. Marsch, in Reviews in Modern Astronomy, edited by G. Klare (Springer- Verlog, Berlin, 1990), p. 
[16] W. H. Matthaeus and M. L. Goldstein, J. Geophys. Res. 87, 6011 (1982). 
[17] E. Marsch and C.-Y. Tu, J. Geophys. Res. 95, 8211 (1990). 

[18] M. K. Vcrma, D. A. Roberts, and M. L. Goldstein et al, J. Geophys. Res. 101, 21619 (1996). 

[19] W. C. Miiller and D. Biskamp, Phys. Rev. Lett. 84, 475 (2000). 

[20] D. Biskamp and W. C. Miiller, Phys. Plasma 7, 4889 (2000). 

[21] S. Sridhar and P. Goldreich, Astrophys. J. 432, 612 (1994). 

[22] P. Goldreich and S. Sridhar, Astrophys. J. 438, 763 (1995). 

[23] J. D. Fournier, U. Frisch, and A. Rose, J. Phys. A 11, 187 (1978). 

[24] W. D. McComb and A. G. Watt, Phys. Rev. A 46, 4797 (1992). 

[25] Y. Zhou and G. Vahala, Phys. Rev. E 47, 2503 (1993). 

[26] A. Pouquet, J. Fluid Mech. 88, 1 (1994). 

[27] A. Ishizawa and Y. Hattori, J. Phy. Soc. Jpn., 67, 4302, (1998). 

[28] A. Ishizawa and Y. Hattori, J. Phy. Soc. Jpn. 67, 441 (1998). 

[29] G. Dar, M. K. Verma, and V. Eswaran, Physica D, 157, 207, 2001. 

[30] M. K. Verma and J. K. Bhattacharjee, Europhys. Lett. 31, 195 (1995). 



14 



Figure Captions 



Fig. 1. Plot of u*(k') (solid) and rj*(k') (dashed) vs. k' for d = 3 and a c = 0, ta = 1- Values at various iterations 
are shown by different curves. 

Fig. 2. Plot of v*(k') and rf{k') vs. k' for d = 2 and cr c = 0,r A = 1. 

Fig. 3. Plot of asymptotic v* (square) and rf (diamond) vs. d for a c = and ta = 1- The fluid ^* (triangle) is 
also plotted for reference. The solid lines are the d^ 1 ^ 2 curves. 
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TABLE I. The values of v* , 77*, v uu * , v ub * , r) bu * , V bb * for various space dimensions d with r A = 1 and a c = 0. 



d 


V* 


V* 


Pr 








r, bb * 


2.1 

2.2 

2.5 

3.0 

4.0 

7.0 

10.0 

50.0 

100 


1.9 
1.2 
1.00 
0.83 
0.62 
0.51 
0.23 
0.14 


0.32 
0.57 
0.69 
0.70 
0.59 
0.50 
0.23 
0.14 


6.0 
2.1 
1.4 
1.2 
1.1 
1.0 
1.0 
1.0 


-0.041 

0.089 

0.20 

0.27 

0.26 

0.23 

0.11 

0.065 


1.96 
1.15 
0.80 
0.56 
0.36 
0.28 
0.12 
0.069 


-0.44 
-0.15 
0.078 
0.21 
0.25 
0.22 
0.11 
0.066 


0.76 
0.72 
0.61 
0.49 
0.34 
0.28 
0.12 
0.069 




TABLE II. 


The values of v* . 




v ub *,ri bu *,ri b '"' for 


various rA when d = 


3 and a c — 0. 




rA 


V* 


V* 


Pr 








r, bb * 


oo 

5000 

100 

5 

2 

1 

0.5 
0.3 


0.38 
0.36 
0.36 
0.47 
0.65 

I. 00 
2.1 

II. 


0.85 
0.85 
0.82 
0.78 
0.69 
0.50 
0.14 


0.42 
0.42 
0.57 
0.83 
1.4 
4.2 
78 


0.38 
0.36 
0.36 
0.32 
0.27 
0.20 
0.11 
0.022 


1.4E-4 
7.3E-3 
0.15 
0.38 
0.80 
2.00 
11.0 


-0.023 
-0.022 
-4.7E-4 
0.031 
0.078 
0.15 
0.082 


0.87 
0.87 
0.82 
0.75 
0.61 
0.35 
0.053 


0.2 

















1G 




FIG. 1. Plot of v*{k') (solid) and if(k') (dashed) vs. k' for d — 3 and a c — Q,ta~ 1. Values at various iterations are shown 
by different curves. 
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FIG. 2. Plot of v*(k') and rf{k') vs. k' for d = 2 and <t c = 0, rA = 1. 
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d 



FIG. 3. Plot of asymptotic v* (square) and (diamond) vs. d for a c — and ta = 1. The fluid j/* (triangle) is also plotted 
for reference. The solid lines are the gT 1 ' 2 curves . 
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